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Abstract 

A+- and A^-hypernuclei are studied in the quark-meson coupling (QMC) model. Com- 
parisons are made with the results for A-hypernuclei studied in the same model previously. 
Although the scalar and vector potentials felt by the A, A+ and Aj, in the corresponding 
hypernuclei multiplet which has the same baryon numbers are quite similar, the wave 
functions obtained, e.g., for lsj/2 state, are very different. The A+ baryon density dis- 
tribution in ^+Pb is much more pushed away from the center than that for the A in 

^ 09 Pb due to the Coulomb force. On the contrary, the Aj, baryon density distributions in 
A;,-hyper nuclei are much larger near the origin than those for the A in the corresponding 
A-hypernuclei due to its heavy mass. It is also found that level spacing for the Aj, single- 
particle energies is much smaller than that for the A and A+. 
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Recently we have made a systematic study of the changes in properties of the heavy 
hadrons which contain a charm or a bottom quark in nuclear matter [Q]. The results suggest 
that the formations of charmed and bottom hypernuclei, which were predicted first in mid 
70's [0, H|, are quite likely. The experimental possibilities were also studied later [|J. In 
addition, we predicted the 5~-nuclear (atomic) bound states, based on analogy with kaonic 
atom H, and a study made for the D- and -D-nuclear bound states || using the quark-meson 
coupling (QMC) model |], |, |. 

The QMC model, which was used there and is used in this study, has been successfully 
applied to many problems associated with nuclear physics and hadronic properties in nuclear 
medium |I| [□]]. For example, the model was applied to the study of J/\l/ dissociation in nuclear 
matter ||12|| , and D and D productions in antiproton-nucleus collisions |JL3[| . Furthermore, 
although only limited studies for heavy mesons (not for heavy baryons) with charm in nuclear 
matter were made by QCD sum rule (for J/ty [|TJ], [TS| and D(D) [0), a study PJ for heavy 
baryons with a charm or a bottom quark based on quarks was made using QMC. In particular, 
recent measurements of polarization transfer performed at MAMI and Jlab |I7]] support the 
medium modification of the proton electromagnetic form factors calculated by the QMC model. 
The final analysis seems to become more in favor of QMC, although still error bars may 
be large to draw a definite conclusion. 

Certainly, the model has shortcomings to be improved eventually. Difficulties to handle it 
will be increased rapidly if we adopt Hartree-Fock approximation even for nuclear matter [|19| . 
and the inclusion of Pauli blocking at the quark level, and the EiV — AiV channel coupling 
have not been implemented yet in a consistent manner with the underlying quark degrees of 
freedom [|l(| . (It should be mentioned that in the case of S-hypernuclei no narrow states have 
been observed. It is unlikely that it will be possible to find such states in the present case |]20|| .) 
Furthermore, an application to double hypernuclei has not been attempted, although recently 
the existence was confirmed ||21|| . Nevertheless, with its simplicity and successful applicability 
achieved so far, we feel some confidence that such a quark-meson coupling model will provide 
us with valuable glimpse into the properties of charmed- and bottom-hypernuclei. 

In this article, we make a quantitative study for the A+- and A^-hypernuclei, by solving 
a system of equations for finite nuclei, embedding a A+ or a A& to the closed-shell nucleus 
in Hartree, mean-field, approximation. Then, the results are compared with those for the A- 
hypernuclei |l(J, which were studied in QMC. It is shown that, although the scalar and vector 
potentials felt by the A, A+ and A& in the corresponding hypernuclei multiplet which has the 
same baryon numbers are quite similar, the wave functions obtained, e.g., for lsi/2 state are 
very different. Namely, the Aj" baryon density distribution in ^°+Pb is much more pushed away 
from the center than that for the A in ^ 09 Pb due to the Coulomb force. On the contrary, the 
Ab baryon density distributions in Afe-hypernuclei are much more central than those for the A 
in the corresponding A-hypernuclei due to its heavy mass. In addition it turns out that the 
level spacing for the A^ single-particle energies is much smaller than that for the A and A+, 
which may imply many interesting new phenomena, which will be discovered in due course by 
experiments. We hope this study opens a new possibility for experiments, related to nuclear 
and hadronic physics, especially for Japan Hadron Facility (JHF). 

We start to consider static, (approximately) spherically symmetric charmed and bottom 
hypernuclei (closed shell plus one heavy baryon configuration) ignoring small non-spherical 
effects due to the embedded heavy baryon. We adopt Hartree, mean-field, approximation. In 



1 



this approximation, pNN tensor coupling gives a spin-orbit force for a nucleon bound in a 
static spherical nucleus, although in Hartree-Fock it can give a central force which contributes 
to the bulk symmetry energy || |§. Furthermore, it gives no contribution for nuclear matter 
since the meson fields are independent of position and time. Thus, we ignore the pNN tensor 
coupling in this study as usually adopted in the Hartree treatment of quantum hadrodynamics 
(QHD) H. 

Using the Born-Oppenheimer approximation, a relativistic Lagrangian density which gives 
the same mean-field equations of motion for a charmed or a bottom hypernucleus, in which 
the quasi-particles moving in single-particle orbits are three-quark clusters with the quantum 
numbers of a charmed baryon, a bottom baryon or a nucleon, when expanded to the same order 
in velocity, is given by QMC R g, p 
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M c (v) ~ ( <7>-(r) + g c p I^b(r) + eQ c A(r) ) 7 o ^c(r), (1) 



where ^(r) (^(r)) and 6(r) are respectively the nucleon (charmed and bottom baryon) and 
the p meson (the time component in the third direction of isospin) fields, while m CT , m w and 
m p are the masses of the a, uj and p meson fields, and g p are the cj-A^ and p-N coupling 
constants which are related to the corresponding (u, d)-quark-u;, g%, and (u, d)-quark-p, g q , 
coupling constants as g u = 3g^ and g p = g q @, . Hereafter we will use notations for the quark 
flavors, q = u,d and Q = s, c, b. 

In an approximation where the a, uj and p fields couple only to the u and d quarks, the 
coupling constants in the charmed and bottom baryons are obtained as g% = (n q /3)g w , and 
g p = g p = g p , with n q being the total number of valence u and d (light) quarks in the baryon 
C. 1% and Qc are the third component of the baryon isospin operator and its electric charge in 
units of the proton charge, e, respectively. The field dependent a-N and a-C coupling strengths 
predicted by the QMC model, g a (cr) and g^(cr), related to the Lagrangian density, Eq. (0), at 
the hadronic level are defined by: 
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M N 
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where {Mc) is the free nucleon (charmed and bottom baryon) mass (masses). Note that 
the dependence of these coupling strengths on the applied scalar field must be calculated 
self-consistently within the quark model || |9|, |10[| . Hence, unlike quantum hadrodynamics 
(QHD) [p2| , even though g% (a) / g a (a) may be 2/3 or 1/3 depending on the number of light 
quarks in the baryon in free space (a = 0)0, this will not necessarily be the case in nuclear 



1 Strictly, this is true only when the bag radii of nucleon and baryon C are exactly the same in the present 
model. See the last line in Eq. (pTl). 



2 



matter. More explicit expressions for g%(cr) and g a (cr) will be given later. From the Lagrangian 
density, Eq. ([!]), a set of equations of motion for the charm or bottom hypernuclear system is 
obtained, 



[i 7 • d - M%{a) - ( gjutf + g P ^b(r) + -(1 + r 3 N )A(f) ) 7 o]^(r) = 0, (4) 
[i 7 • d - Mq(g) - ( gfair) + g p I^b(r) + eQ c A(r) ) 7 o]^c(r) = 0, (5) 
(-V 2 r + m 2 a )a(^ = - ? K > \ Ps (r) - [ ^ > f(r), 

= g„C N {<T)p.{r) + g c a C c {(j)p c s {?)i (6) 

(-V r 2 + m 2 , ) w (r) = &,pb(tO + <&g(r), (7) 

("V r 2 + m>(r) = |p 3 (r) + gf I^p c B (f), (8) 

(-V 2 )A(r) = e Pp (r) + eQ c p C B (f), (9) 

where, p s {f) (Ps(r)), Pb{?) {Pb{^))i Pzif) an d Ppif) are the scalar, baryon, third component of 
isovector, and proton densities at the position r in the charmed or bottom hypernuclei || ||, [[(J 



On the right hand side of Eq. ©, -[^%&] = ^Cjv(a) and -[^^] = g°C c {jj), where 
g CT = g CT (cr = 0) and g% = g„(<J = 0), are a new, and characteristic feature of QMC beyond 



QHD p2| |2J| . The effective mass for the charmed or bottom baryon C is defined by, 

-n,<£ £ dyi> q {y)i> q {y) = -n q glS c {a) = -— [g^(a)a] , (10) 
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with the MIT bag model quantities [7L 0, P, 10], 
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CcO) = Sc(<t)/S c {0), 9° = n q g*S c (0) = -^g„S c (0)/SW(0) = -^^r c/iV . (11) 



Quantities for the nucleon are similarly obtained by replacing the indices, C — > iV '. Here, 
S, # g ,Q, and m 9i Q are the parameters for the sum of the cm. and gluon fluctuation effects, bag 
pressure, lowest eigenvalues for the quarks, q or Q, respectively, and the corresponding current 
quark masses. Zn and B [zc) are fixed by fitting the nucleon (charmed or bottom baryon) mass 
in free space. Concerning the sign of m* in (hyper)nucleus, it reflects nothing but the strength 
of the attractive scalar potential, and thus naive interpretation of the mass for a (physical) 
particle, which is positive, should not be applied. 

The bag radii in- medium, R* NC , are obtained by the equilibrium condition 
dM^ c (a) /dR]y : c\R N C =R* N c = 0- The bag parameters calculated and chosen for the present 
study in free space are, (z^, z^, z A +, z\ b ) = (3.295,3.131,1.766,-0.643), (Rn, R\, R A +, Rx b ) 
= (0.800,0.806,0.846,0.930) fm, B 1 ^ = 170 MeV, (m u4 ,m s ,m c ,m b ) = (5,250,1300,4200) 
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MeV. The parameters associated with the u, d and s quarks are the same as in our previous 
investigations || [K|. At the hadron level, the entire information on the quark dynamics is 
condensed in Cn,c(°') of Eq. (||). The parameters at the hadron level, which are already 
fixed by the study of infinite nuclear matter and finite nuclei ||, are as follows: = 783 
MeV, m p = 770 MeV, m a = 418 MeV, e 2 /4n = 1/137.036, g^/in = 3.12, gl/in = 5.31 and 
^/4vr = 6.93. 

We briefly discuss about the spin-orbit force in QMC ||. The origin of the spin orbit 
force for a composite nucleon moving through scalar and vector fields which vary with position 
was explained in detail in Ref. || - c.f. sect. 3.2. The situation for the A and also for other 
hyperons are discussed in detail in Ref. [|lOfl . 

In order to include the spin-orbit potential (approximately) correctly, e.g., for the A+, we 
add perturbatively the correction due to the vector potential, — 2M J ij-g^ c uj(r) \ I ■ s, to 

the single-particle energies obtained with the Dirac equation, in the same way as that added 
in Ref. ||10|| . This may correspond to a correct spin-orbit force which is calculated by the 
underlying quark model R CT : 

V &W- *= -^^y (ll M h(?) + ^(Ol) r- s, (12) 
since the Dirac equation at the hadronic level solved in usual QHD-type models leads to: 

V&W- s = (liK+irl - f. s, (13) 

which has the opposite sign for the vector potential, g^u{f). The correction to the spin-orbit 
force, which appears naturally in the QMC model, may also be modeled at the hadronic level 



of the Dirac equation by adding a tensor interaction, motivated by the quark model |j25 , 2£|. 
Here, we should make a comment that, as was discussed by Dover and Gal p?| in detail, one 
boson exchange model with underlying (approximate) SU(3) symmetry in strong interaction, 
also leads to the weaker spin-orbit forces for the (strange) hyperon-nucleon (YN) than that for 
the nucleon-nucleon (NN). 

However, in practice, because of its heavy mass (M*+), contribution to the single-particle 
energies from the spin-orbit potential with or without including the correction term, turned out 
to be even smaller than that for the A-hypernuclei, and further smaller for the Ab-hypernuclei. 
Contribution from the spin-orbit potential with the correction term is typically of order 0.01 
MeV, and even for the largest case is = 0.1 MeV. This can be understood when one considers 
the limit, — ► oo in Eq. (|T2|), where the quantity inside the square brackets varies smoothly 
from an order of hundred MeV to zero near the surface of the hypernucleus, and the derivative 
with respect to r is finite. (See also Figs. |2] and |3|.) 

Now we discuss the results. First, we show in Fig. [I] the total baryon density distributions 
in ^Ca and ^ 09 Pb (j = A, A+,A&), for lsi/ 2 configuration in each hypernucleus. Note that 
because of the self-consistency, the total baryon density distributions are dependent on the 
configurations of the embedded particles. The total baryon density distributions are quite 
similar for the A-, A+- and Afe-hypernuclei multiplet which has the same baryon numbers, A, 
since the effect of A, A+ and A& is = 1/A for each hypernucleus. Nevertheless, one notices that 
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the A^-hypernuclei density near the center is slightly higher than the corresponding A- and A+- 
hypernuclei. This is because the is heavy and localized nearer the center, and contributes 
to the total baryon density there. The baryon (probability) density distributions for the A, A+ 
and Af, in corresponding hypernuclei will be shown later. 

Next, in Figs. |2] and ^|, we show the scalar and vector potentials felt by the A, A+ and 
A b for lsi/2 state in f'Ga and ^ 09 Pb (j = A, A+, A b ), and the corresponding probability density 
distributions in Fig. 4. In Figs. 2 and 3 "Pauli" stands for the effective, repulsive, potential 
representing the Pauli blocking at the quark level plus the T, C)b N — A Cjb N channel coupling, 
introduced at the baryon level phenomenologically |I(J. For the A+, the Coulomb potentials 
are also shown. As for the case of the nuclear matter |IJ, the scalar and vector potentials 
felt by these particles in hypernuclei multiplet which has the same baryon numbers are also 
quite similar. Thus, as far as the total baryon density distributions and the scalar and vector 
potentials are concerned, A-, A+- and A(,-hypernuclei show quite similar features within the 
multiplet. However, as shown in Fig. 4, the wave functions obtained for lsi/2 state are very 
different. The A+ baryon density distribution in ^°+Pb is much more pushed away from the 
center than that for the A in ^ 09 Pb due to the Coulomb force. On the contrary, the A b baryon 
density distributions in A^-hypernuclei are much larger near the origin than those for the A in 
the corresponding A-hypernuclei due to its heavy mass. 

Having obtained reasonable ideas about the potentials felt by A, A+ and A&, we show the 
calculated single-particle energies in Tables [I] and 0. Results for the A-hypernuclei are from 
Ref. ||10|| . In these calculations, effective Pauli blocking, the effect of the Y; Cib N — A C}b N channel 



coupling, and correction to the spin-orbit force based on the underlying quark structure are 
included in the same way as adopted in Ref. [|10j- (However, recall that the negligibly small 
contribution from the correction terms for the spin-orbit force, and also contributions from the 
spin-orbit force itself.) Note that since the mass difference of the A+ and S c is larger than 
that of the A and E, and it is probably also true for the A b and we expect the effect of the 
channel coupling for the charmed and bottom hypernuclei is smaller than those for the strange 
hypernuclei, although the same parameters were used in the present calculation. In addition, we 
searched for the single-particle states up to the same highest state as that of the core neutrons 
in each hypernucleus, since the deeper levels are usually easier to observe in experiment. 

In Tables [1] and |2|, it is clear that the A+ single-particle energy levels are higher than 
the corresponding levels for the A and A;,. This is a consequence of the Coulomb force. This 
feature becomes stronger as the proton number in the core nucleus increases. 

Second, the level spacing for the A b single-particle energies is much smaller than that 
for the A and A+. This may be ascribed to its heavy mass (or M b ). In the Dirac equation 
for the Ab, the mass term dominates more than that of the term proportional to Dirac's k, 
which classifies the states, or single-particle wave functions. (See Refs. [[|, |l(J for detail.) This 
small level spacing would make it very difficult to distinguish the states in experiment, or to 
achieve such high resolution. On the other hand, this may imply also many new phenomena. It 
will have a large probability to trap a A^ among one of those many states, especially in heavy 
nucleus such as lead (Pb). What are the consequences ? May be the A^ weak decay happens 
inside a heavy nucleus with a very low probability ? Does it emit many photons when the A b 
gradually makes transitions from a deeper state to a shallower state ? All these questions raise 
a flood of speculations. 

To summarize, we have made a quantitative study for A+- and A;,-hypernuclei in a quark- 
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meson coupling model. We have solved a system of equations in a self-consistent approach for 
several finite nuclei with closed shell plus a hyperon (A+ or A&), embedding a A+ or A b in the 
nucleus. Results are compared with those for the A-hypernuclei. It is shown that, although the 
scalar and vector potentials felt by the A, A+- and A b are quite similar in corresponding hyper- 
nuclei multiplet which has the same baryon numbers, the single-particle wave functions, and 
single-particle energy level spacings are quite different. For the A+-hypernuclei, the Coulomb 
force plays a crucial role, and so does the heavy A b mass for the A 6 -hypernuclei. It should be 
emphasized that we have used the values for the coupling constants of er (or a-field dependent 
strength), uj and p to the A, A+ and A b determined automatically based on the underlying quark 
model, as for the nucleon and other baryons. (Recall that the values for the vector u fields to 
any baryons can be obtained by the number of light quarks in a baryon, but those for the a are 
different as shown in Eqs. (|i0|) and fllTD -) Phenomenology would determine ultimately if the 
coupling constants (strengths) determined by the underlying quark model actually work for A+ 
and A b or not. Although implications of the present results can be speculated a great deal, we 
would like to emphasize that, what we showed is that the A+- and Ab-hypernuclei would exist 
in realistic experimental conditions. Experiments at facilities like JHF would provide further 
inputs to gain a better understanding of the interaction of A+ and A b with the nuclear matter. 
Additional studies are needed to investigate the semi-leptonic weak decay of A+ and A b hyper- 
ons. The role of Pauli blocking and density in influencing the decay rates as compared to the 
free hyperons would be highly useful. Such study can have an impact on the hadronization of 
the quark-gluon plasma and the transport of hadrons in nuclear matter of high density. Will 
the high density lead to a slower decay and that a higher probability to survive its passage 
through the material ? At present the study of the presence of A+ and A b in finite nuclei is 
its infancy. Careful investigations, both theoretical and experimental, would lead to a much 
better understanding of the role of heavy quarks in finite nuclei. 
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Table 1: Single-particle energies (in MeV) for ) 7 0, f Ca and f Ca (j = A,A+,A 6 ). Sin- 
gle-particle energy levels are calculated up to the same highest states as that of the core neu- 
trons. Results for the hypernuclei are taken from Ref . [ IU| . Experimental data for A-hypernuclei 



are taken from Ref. 
by the experiments. 



where spin-orbit splittings for A-hypernuclei are not well determined 





A U 

(Exp.) 


17 D 

A u 


17 D 


17 D 

A b U 


fCa 
(Exp.) 


^ Ccl 


A+ Ca 




Ccl 




49 p 


lSl/2 


-12.5 


-14.1 


-12.8 


-19.6 


-20.0 


-19.5 


-12.8 


-23.0 


-21.0 


-14.3 


-24.4 


1^3/2 


-2.5 


-5.1 


-7.3 


-16.5 


-12.0 


-12.3 


-9.2 


-20.9 


-13.9 


-10.6 


-22.2 




(IPs/a) 


-5.0 


-7.3 


-16.5 


(IP3/2) 


-12.3 


-9.1 


-20.9 


-13.8 


-10.6 


-22.2 


14/2 












-4.7 


-4.8 


-18.4 


-6.5 


-6.5 


-19.5 


2si/ 2 












-3.5 


-3.4 


-17.4 


-5.4 


-5.3 


-18.8 


14/2 












-4.6 


-4.8 


-18.4 


-6.4 


-6.4 


-19.5 


1/7/2 




















-2.0 


-16.8 
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Table 2: Single-particle energies (in MeV) for ^Zr and | 08 Pb (j = A, A+,Af>). Experimental 
data are taken from Ref. ]29]. See caption of Table |l]for other explanations. 





A 9 Yb 
(Exp.) 


91 r7 


Q1 rj 

A+ Zr 


91 r7 


2(18 T-11 

A U8 Pb 
(Exp.) 


209 ni 


209 ni 

f+Pb 

**-c 


209 ni 
209p b 




-22.5 


-23.9 


-10.8 


-25.7 


-27.0 


-27.0 


-5.2 


-27.4 


Ipa/a 


-16.0 


-18.4 


-8.7 


-24.2 


-22.0 


-23.4 


-4.1 


-26.6 


lPl/2 


(IP3/2) 


-18.4 


-8.7 


-24.2 


O3/2) 


-23.4 


-4.0 


-26.6 


14/2 


-9.0 


-12.3 


-5.8 


-22.4 


-17.0 


-19.1 


-2.4 


-25.4 


2si/ 2 





-10.8 


-3.9 


-21.6 





-17.6 





-24.7 


14/2 


(W5/2) 


-12.3 


-5.8 


-22.4 


(W5/2) 


-19.1 


-2.4 


-25.4 


1/7/2 


-2.0 


-5.9 


-2.4 


-20.4 


-12.0 


-14.4 


— 


-24.1 


2p 3 /2 


— 


-4.2 


— 


-19.5 


— 


-12.4 


— 


-23.2 


1/5/2 


(I/7/2) 


-5.8 


-2.4 


-20.4 


(I/7/2) 


-14.3 


— 


-24.1 


2j?i/a 




-4.1 




-19.5 




-12.4 




-23.2 


109/2 








-18.1 


-7.0 


-9.3 




-22.6 


lflY/2 










(1^9/2) 


-9.2 


— 


-22.6 


1^11/2 












-3.9 




-21.0 


^"5/2 












-7.0 




-21.7 


24/2 












-7.0 




-21.7 


l/ig/2 












-3.8 




-21.0 


3S1/2 












-6.1 




-21.3 


2/7/2 












-1.7 




-20.1 


3^3/2 












-1.0 




-19.6 


2/5/2 












-1.7 




-20.1 


3^1/2 












-1.0 




-19.6 


H13/2 
















-19.3 
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Figure 1: Total baryon density distributions in ^Ca and | 09 Pb (j = A, A+,A & ), for lsi/ 2 
configuration for the A, A+ and A b . 
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Figure 2: Potential strengths for lsx/ 2 state felt by the A,A+ and A b in ^Ca (j = A,A+,A fe ). 
"Pauli" stands for the effective, repulsive, potential representing the Pauli blocking at the quark 
level plus the T> C ^N — A C ^N channel coupling, introduced at the baryon level phenomenologi- 
cally 0. 
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Figure 3: Potential strengths for lsi/ 2 state felt by the A, A+ and A b in ^ 09 Pb (j = A, A+, A fe ). 
See also caption of Fig. ||. 
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Figure 4: A,A+ and A b baryon (probability) density distributions for lsi/ 2 state in ^Ca and 
f»Fb(j=A,A+ A 6 ). 
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